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ABSTRACT 



Q-i! 



Context. HD 110432 was classified as a "y Cas X-ray analog" since it has similar peculiar X-ray and optical characteristics, i.e. a 
hard-thermal X-ray variable emission and an optical spectrum affected by an extensive disk. Lopes de Oliveira et al. (2007) suggest 
that it might be a Be star harboring an accreting white dwarf or that the X-rays may come from an interaction between the surface of 
' , ' the star and its disk. 

Aims. To investigate the disk around this Be star we used the VLTI/ AMBER instrument, which combines high spectral (R= 12000) 
' and high spatial (# ra j n =4 mas) resolutions. 

Methods. We constrain the geometry and kinematics of its circumstellar disk from the highest spatial resolution ever achieved on this 
star. 

Results. We obtain a disk extension in the Bry line of 10.2 D* and 7.8 D* in the He I line at 2.05 yum assuming a Gaussian disk model. 
The disk is clearly following a Keplerian rotation. We obtained an inclination angle of 55°, and the star is a nearly critical rotator 
with V rot /V c =1.00±0.2. This inclination is greater than the value found for yCas (about 42°, Stee et al. 2012), and is consistent with 
. the inference from optical Fe II emission profiles by Smith & Balona (2006) that the inclination should be more than the y Cas value. 

In the near-IR continuum, the disk of HD 110432 is 3 times larger than y Cas's disk. We have no direct evidence of a companion 
around HD 110432, but it seems that we have a clear signature for disk inhomogeneities as detected for £ Tau. This asymmetrical 
disk detection may be interpreted within the one-armed oscillation viscous disk framework. Another finding is that the disk size in 
the near-IR is similar to other Be stars with different spectral types and thus may be independent of the stellar parameters, as found 
^ for classical Be stars. 

t> 

■ Key words. Techniques: high angular resolution - Techniques: interferometric - Stars: emission-line, Be - Stars: winds, outflows 

OO ' Stars: individual (y Cas) - Stars: circumstellar matter 

^ . 

C^) 1 . Introduction for a blue straggler. In addition, the star is situated near and be- 

' Tm „ „ , , TTT „, ^ , , yond the Southern Coalsack dark nebula. However, only the UV 

52«1 °f [ ? ?! = InS?' B0 - 5 " 1 m " IVe; 6t aL "eddening and extinction are appreciably affected by the pres- 

(1986), Smith & Balona (2006)) is an interesting object since ence of ^ nebula ^ ^ f ound 

it exhibits variable hard-thermal X-rays, which is very different 

from other "classical" massive stars that are emitting soft X-rays The X-ray properties of HD 110432 have been described 

or from non-thermal emission of all known Be/X ray binaries, well by a series of papers by Lopes de Oliveira et al. (2007), 

, The other well known Be star sharing these properties is y Cas. Torrejon et al. (2012), and Smith et al. (2012). Although the X- 

Smith & Balona (2006) propose that HD 110432 should be the ray flaring properties are very comparable, the last two papers 

first new member of a select new class of "y Cas X-ray analogs." demonstrate that the dominant hot thermal plasma component 

These analogs exhibit the same peculiar X-ray and optical char- associated with HD 1 10432 emission is even hotter than the anal- 

acteristics, i.e. Balmer and Fe II double-peaked emission lines ogous component of y Cas. 

with a strong IR excess due to free-free and free-bound emis- . . _ , _ . , , 

„ , , lrioo o c inoi\ /-i if The origin of this X-Ray emission of these y Cas variables 
sion (Stellebak 1982; Meyer & Savage 1981), as y Cas itself. f . . . . ^_ , ' 
,„,,„,„. J . , f- tt is still a subject of debate. Currently there are two competing 
HD1 10432 is not yet known to be in a binary system. However, it _,, J _ . , , : , _ •■ , 
u u t .L i » KTn^/tnn/r- » • o ./ scenarios. The first one is based on single Be stars with unusual 
may be a member of the cluster NGC 4609 (Feinstein & Marraco . . . , . , , . 
, mm T .. ,-niv/r a ' *■ ' a- a i strong magnetic activity and rapid X-ray correlations with op- 
1979). In this case its age would be 60 Myr, and it is a candidate . , & , * , . . D _ A TT o 
J tical and UV activity (see Smith & Robinson (2003), Henry & 

Send offprint requests to: Philippe.Stee@oca.eu Smith (2012), Smith et al. (2012)). The second scenario proposes 

* Based on observations made with VLTI ESO telescopes at La a binary system with a Be star and an accreting degenerate com- 

Silla Paranal Observatory under GTO programme IDs 0.84.C-0062(A), panion, such as a white dwarf in accordance with the evolution- 

0.84.C-0062(B), 0.84.C-0062(C), 0.84.C-0062(D) ary models of massive binary systems. 
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To better constrain the physics and the mechanism respon- 
sible for this X-ray emission, as well as the link to its circum- 
stellar environment, we observed this Be star with the Very 
Large Telescope Interferometer (VLTI) combining 1.8m tele- 
scopes on baselines up to 128m, with the AMBER focal instru- 
ment that combines high spectral (R= 12000) and high spatial 
(0min— 4mas) resolutions (Petrov et al. 2007; Robbe-Dubois et 
al. 2007) mainly focused on the Bry and He I 2.05 //m emission 
lines. 

In this paper, we would like to address the following points. 

• If there is a strong coupling between a putative magnetic 
field and the circumstellar disk that may lead to part of the 
observed X-ray emission, is the disk forced to rotate as a 
solid body by the magnetic field? What is the kinematics 
within HD1 10432's circumstellar environment? 

• What is the disk extension in the Bry and He 1 2.05 yum emis- 
sion lines? Is the circumstellar disk dense and/or large as 
suggested from the strong and quasi-symmetrical profile of 
the Ha line (EW ~ 60 A) and the detection of several metal- 
lic lines in emission by Lopes de Oliveira et al. (2007) and 
Torrejon et al. (2012)? 

• What is the central Be star's rotational velocity? Is the stellar 
rotation close to critical as already found by Meilland et al. 
(2012) for eight other Be stars? 

• Since the binarity of HD 1 10432 may figure in the produc- 
tion of the hard-thermal X-ray emission, either as a source of 
X-ray emission or previous transfer of angular momentum to 
the Be star, do we have any evidence of a companion in our 
interferometric data? 

The paper is organized as follows. In Section 2 we summarize 
the stellar parameters of HD 110432. In Section 3 we present 
the VLTI/ AMBER observations and the data reduction process. 
A first analysis using an axi-symmetrical geometrical models is 
presented in Sects. 4 and 5, and a more advanced modeling of the 
differential data using a kinematic model is discussed in Section 
6. Finally, we summarize our results in Section 7. 

2. The stellar parameters and distance 

To constrain the stellar contribution to the total flux we need 
to know the effective temperature T e ff, the radius and the 
gravity log g e g of the central star. However, these stellar pa- 
rameters for rapidly rotating Be stars are quite uncertain. Taking 
gravitational darkening effects into account, Fremat et al. (2005) 
showed that the apparent T e g is typically thought to be ~ 5 % to 
15% colder than the values of the non rotating counterparts. 
In the case of HD 110432, they found an apparent T e $ of 
20324+344K and T eS of 24070+603 for its parent non rotat- 
ing counterpart (pnrc). Using the same method, Zorec et al. 
(2005) derived an average T e g on the stellar surface of 225 10 K. 
Knowing these uncertainties on the definition of the effective 
temperature, we adopted T e g =22000+2000K in this paper. 
On the other hand, Fremat et al. find an apparent g e ff of 
3.638+0.042 and 3.950+0.074 for the pnrc. As this parameter 
does not strongly affect the overall SED and the determination of 
the stellar fluxes in the infrared, we simply assume g e ff =3.9+0. 1 . 
The stellar radius is also strongly affected by the stellar rotation. 
Zorec et al. (2005) then estimate the stellar mass as M/M Q = 9.6 
and the stellar age f/fjyis = 0.61 (fMS is the time the rotating 
star spend in the main sequence), which mean that the corre- 
sponding 'spherical' radius of the parent non rotating object is 
R- k =5.3R , which leads to the critical radius at the equator of 



Parameter 


T c ff 


log geff 


R* 


d 


(unit) 


(K) 




(Ro) 


(pc) 


Value 


22000 


3.9 


6.5 


373 


Uncertainty 


+2000 


±0.1 


±1.2 


±41 



Table 1. Stellar parameters and distance of HD 1 10432 adopted 
in this paper. 



TJRq. We note that the mean apparent radius will depend on 
the object inclination angle. However, here we decide to adopt 
a value of /?+=6.5±1.2R Q . Using the distance derived from van 
Leeuven (2007) Hipparcos parallaxes, i.e. d=373+pc, we infer a 
photosphere angular size of 0.16 mas. The parameters are sum- 
marized in Table Q] 

Using this set of stellar parameters and photometric measure- 
ments from the SIMBAD Q database, we tried to determine 
the relative flux of the circumstellar environment within the 
VLTI/ AMBER spectral domain, i.e. the H and K bands. To de- 
termine the stellar flux in these bands we used Kumcz models 
reddened using the law of extinction from Cardelli et al. (1989). 
However, we found that the influence of the reddening in the 
near-infrared fluxes is negligible compared to the uncertainties 
on the stellar parameters and distance. 

Finally we found that the emission in the H and K bands are 
both dominated by the circumstellar environment with a relative 
circumstellar flux (Fdi s k) of 71+7% and 79±5%, respectively. 



3. Observations and data reduction process 

HD 110432, was observed in January, February and March 
2010 with the A0-K0-G1 and D0-H0-G1 VLTI/ AMBER triplets 
of the 1.8m Auxiliary Telescopes array, ranging from 51 to 
124m baselines (See Fig. [TJ. We used both low resolution 
(LR=30) and high resolution (HR= 12000) measurements in the 
H ( 1 .54- 1 .87//m) and K ( 1 .94-2.37/mi) bands with the help of the 
FINITO fringe tracker. With an apparent diameter of 0.83 ± 0.01 
mas in the CHARM2 catalog from Richichi et al. (2005), HD 
113752 was used to calibrate the visibility measurements (see 
Table 2). More details about the AMBER instrument can be 
found in Petrov et al. (2007). 

Data were reduced using the VLTI/AMBER data reduction 
software, i.e., amdlib v3.Q.3bl (see Tatulli et al. 2007 and 
Chelli et al. 2009 for detailed information on the AMBER data 
reduction). We selected individual exposures with the standard 
selection criteria (Millour et al. 2007). We rejected the 80% of 
frames with the lowest S/N. For observations in LR mode we 
also rejected the frames with a piston larger than 10/im, as well 
as frames with a flux ratio between the beams higher than three. 

The interferometric observables (visibility, differential 
phase, and closure phase) were then averaged and calibrated. For 
this last step we used scripts described in Millour et al. (2007) 
that are now part of the standard amdlib package. The calibra- 
tion process includes an estimation of the calibrators' size and 
their uncertainties from various catalogs, a determination of the 
transfer functions and their evolution during the whole night, and 
a computation of the calibrated visibilities and phases. The final 
errors on the measurements include uncertainties on the calibra- 
tors' diameter, the atmosphere transfer function fluctuations, and 
intrinsic errors on the measurements (Fig. |2}. 



1 available at http : / /simbad . u- strasbg . f r/simbad/ 
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Observation Projected baseline Mode DIT Seeing Calibrators 



Date 


Time 


Triplet 


T (m\ 


PA i°\ 






\ ) 




24/01/2010 


07 


32 


A0-K0-G1 


78.7/81.1/127.9 


-163.2/-89.5/-125.7 


T T T~l T r | — i r\ in 

HR-K-F-2.17 


12 


0.69 


113752 


26/01/2010 


07 


33 


A0-K0-G1 


78.5/81.8/127.8 


ly'l n / fin A 1 1 T> ^"7 

-161.7/-87.4/-123.7 


T TT1 T Z | — * r-\ in 

HR-K-F-2.17 


12 


0.60 


113752 


10/02/2010 


07 


33 


A0-K0-G1 


76.1/86.0/126.2 


-151.5/-73.6/-109.7 


HR-K-F-2.17 


12 


0.66 


113752 


10/02/2010 


08 


29 


A0-K0-G1 


73.2/88.3/123.5 


1/11 1 / ZT^ J / AO 1 

-143.1/-62.4/-98.2 


HR-K-F-2.06 


12 


0.84 


1 13752 


1 //Uj/2010 


07 


30 


UU-rtU-Ol 


68.2/55.5/50.5 


-31.8/95.9/24.1 


T 1 > III' I 


U.U5 


U.J6 


1 1 5 Ijl 


17/03/2010 


07 


41 


D0-H0-G1 


68.4/58.0/53.1 


-29.8/101.5/25.5 


LR-HK-F 


0.05 


0.50 


113752 


17/03/2010 


08 


15 


D0-H0-G1 


69.0/56.1/51.5 


-23.1/109.6/30.0 


LR-HK-F 


0.05 


0.57 


113752 


19/03/2010 


05 


23 


A0-K0-G1 


75.6/86.5/125.8 


-149.8/-71.4/-107.4 


LR-HK-F 


0.05 


0.41 


113752 


19/03/2010 


06 


00 


A0-K0-G1 


73.6/88.1/123.9 


-144.0/-63.7/-99.5 


LR-HK-F 


0.05 


0.37 


113752 


19/03/2010 


06 


34 


A0-K0-G1 


71.3/89.2/121.5 


-138.7/-56.4/-92.0 


LR-HK-F 


0.05 


1.30 


113752 


19/03/2010 


07 


08 


A0-K0-G1 


68.6/89.9/118.5 


-133.4/-49.2/-84.4 


LR-HK-F 


0.05 


0.90 


113752 



Table 2. Observations log of HD 1 10432 
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Fig. 1. (u,v) plane coverage for our VLTI/AMBER observations 
of HD 1 10432. LR observations are plotted as red squares, HR 
data as blue circles for Bry, and green ones for the 2.06 //m He I 
line. 



4. Modeling the H and K band continuum 

4.7. The LITpro software 

To model the continuum visibility modulus and closure phases 
obtained from our VLTI/AMBER observations we used the 
LITprcfl model-fitting software for optical/infrared interfer- 
ometric data developed by the Jean-Marie Mariotti Center 
(JMMC) to analyze our data (Tallon et al. 2008). It is based on 
the Levenberg-Marquardt algorithm, which allows a fit to con- 
verge to the closest^ 2 local minimum from a set of initial values 
of the model parameters. It also includes tools for facilitating 
the search for the global minimum. LITpro calculates an error 
on the fitted parameters based on the x 2 value at the minimum. 
It uses data error estimates based on the OI FITS format, which 
does not include error correlation estimates. Therefore, in some 
cases, LITpro can provide underestimated errors on the parame- 
ters. 



2 LITpro software available at http://www.jmmc.fr/litpro 



4.2. Extension of HD1 10432 in the continuum 



Even if the spatial resolution of the VLTI/AMBER instrument is 
about 4 mas, it is not necessary to fully resolve the target, i.e. to 
reach the first zero of the visibility function to estimate its size. 
It is always possible to use a partial resolution of the interferom- 
eter to estimate the size of the target with a model that fits the 
data and thus determines the angular size of targets smaller than 
4 mas. See for instance the review of the interferometric observa- 
tions of rapidly rotating stars by van Belle (2012). Thus, to esti- 
mate the disk extension in the near infrared, we model the central 
star as a uniform disk + an extended elliptical Gaussian distribu- 
tion. This model then has four free parameters: the FWHM of the 
major axis of the Gaussian distribution (6>diskX the relative envi- 
ronment continuum flux (Fdi s k), the flattening ratio of the major 
to the minor Gaussian distribution axis (f), and the position angle 
of the major axis of the disk measured eastward from the north 
(P.A). 

Since we have VLTI/AMBER data in two different spectral do- 
mains, namely the H band between 1 .54- 1.87 //m and the K band 
within 1.94-2.37 //m, we tried to constrain both data sets sepa- 
rately. HD 110432 is a variable Be star, and the data used to 
estimate the circumstellar contribution to the total near-IR con- 
tinuum flux is not contemporaneous to our interferometric mea- 
surements. Nevertheless, using the values of Fdisk obtained in 
section 2, i.e. 71+7% and 79+5% in the H and K bands respec- 
tively, we can infer the disk dimension in the continuum. 

The best LITpro models and the corresponding extensions are 
given in Table [3] We obtain a major Gaussian distribution axis 
of 0.92 mas, a flattening ratio of 1.56 in the H band that, assum- 
ing a flat disk, corresponds to 50 ± 5° with a reduced^- 2 of 4.4 
and a RA. of 22+ 5°. In the K band the disk appears to have sim- 
ilar properties: an angular size of 0.95 mas, a flattening of 1 .62 
corresponding to an inclination inclination angle of 52 ± 5° and 
a compatible RA of 17 ± 5°. These extensions are similar to the 
disk size in the near-IR continuum, i.e. 0.82+0.08 mas, obtained 
for y Cas by Stee et al. (2012). Nevertheless, for y Cas it corre- 
sponds to 1.9 D+ and thus HD 1 10432 exhibits a disk size three 
times larger in the near-IR continuum than y Cas. An intensity 
map in the K band is plotted in Fig. [3] The RA. of ~ 20 ° is nei- 
ther perpendicular nor parallel to the polarization measurement 
of 81° by Yudin (2001). The error bars on these disk measure- 
ments are mainly dominated by the uncertainties on the disk flux 
contribution in the H and K bands. 
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H-Band K-Band 




0.2 0.3 
B/\ (cycles/mas) 



0.2 0.3 
B/A (cycles/mas) 



Fig. 2. Left. H band data from the VLTI/ AMBER instrument for different baselines (blue, red and green lines) with a LITpro model 
over-ploted along the minor (dotted line) and major (plain line) disk axis. Right. Same kind of data but for the K band. 
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Fig. 4. HD 1 10432 differential visibility and phase within the Bry line from our four VLTI/ AMBER HR measurements (red line). 
Each row corresponds to one VLTI/AMBER measurement (3 different baselines). The visibility and phase of the best-fit kinematics 
model is overplotted in green. 



5. Modeling the disk in Bry and Hel emission lines 

5.1. The kinematic model 

To quantitatively constrain the velocity fields in the circum- 
stellar environment of HD 110432, we developed a simple two- 
dimensional kinematic model of a rotating and/or expanding 
equatorial disk. This model has already been used to success- 
fully model classical Be stars (see Delaa et al. 2011; Meilland 
et al. 2011; Meilland et al. 2012) and is described in detail in 
Delaa et al. (2011). The model geometry is completely ad-hoc: 
the star is modeled as a uniform disk and the envelope emission 
in the emission line as an elliptical Gaussian distribution with a 



given FWHM and a flattening due to a projection effect of the 
geometrically thin equatorial disk, i.e., / = l/cos{i), where i is 
the star-disk system's inclination angle. 

The emission maps are then combined with a two- 
dimensional projected velocity map of a geometrically thin ex- 
panding and/or rotating equatorial disk. For each spectral chan- 
nel in the line, an iso-velocity map projected along the line of 
sight is then calculated and multiplied by the whole emission 
map in the line. 

The model parameters can be classified into three categories: 
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Fig. 5. HD 1 10432 differential visibility and phase within the He I line from our four VLTI/ AMBER HR measurements (red line). 
The differential visibility is first plotted with the corresponding phase for 3 different baselines with B=73.20m, PA=-143.06°; 
B=88.32m, PA=-62.36°; and B=123.48m, PA=-98.16°. The visibility and phase of the best-fit kinematics model is overplotted in 
green. 




a (mas, increasing towards East) 

Fig. 3. LITpro model of HD 1 10432 Cas disk in the K band as- 
suming a uniform central star with a surrounding Gaussian disk 

1. The stellar parameters: stellar radius (R+), distance (d), in- 
clination angle (i), and disk major-axis position angle (PA). 

2. The kinematic parameters: rotational velocity (V r ot) as mea- 
sured from the stellar surface, expansion velocity at the pho- 
tosphere (Vo), terminal velocity (V M ), and exponents of the 
expansion (y) and rotation (/?) velocity laws. 

3. The disk emission line parameters: disk FWHM in the line 
(a\) and line equivalent width (EW). 

Table 3. Parameters obtained with LITpro for HD 110432 as- 
suming a 0.31 mas central star (see section 2). 



Parameters 


Hband 


Kband 




1.54- 1.87pm 


1.94-2.37pm 


fdisk (mas) 


0.92 ±0.05 


0.95 ±0.07 


#disk (D*) 


5.7 


5.9 


f 


1.56 +0.05 


1.62 ±0.03 


PA. 


22 +5° 


17 ±5° 


Inclination (i) 


50° 


52° 


F disk (%) 


71 +7 


79 ±5 


xi 


4.4 


3.4 



Notes, ttasi; is the FWHM of the major axis of the Gaussian distribution, 
F disk the relative environment continuum flux, f the flattening ratio of 
the major to the minor Gaussian distribution axis, and PA. the position 
angle of the major axis of the disk measured eastward from the north. 
The estimated inclination angle (i) is given assuming a very thin disk. 



The star distance is taken from van Leeuwen (2007), and R+ is 
derived from the fit of the SED. The nine other parameters are 
free. 

If the disk is directly connected to the stellar surface, the 
rotational velocity (V lot ) should be equal to the stellar rotational 
velocity. However, in some cases, V rot may exceed the stellar 
velocity if the star is not critically rotating and some additional 
momentum is transferred to the circumstellar matter. Finally, we 
considered in our modeling that V rot is a free parameter with a 
higher maximum value equal to the critical velocity (V c ). 

5.2. Model fitting and results 

We have computed several hundred models to constrain the pa- 
rameters, determined the uncertainties, and tried to detect any 
degeneracy or linked parameters. Owing to the large number 
of free-parameters, an automatic model-fitting method would 
have resulted in computing millions of models. Moreover, we 
clearly know each parameter's effect on the visibility and phase 
variations (see Meilland et al. 2012). Consequently, we de- 
cided to perform the fit manually. We could exclude models 
with significant expansion velocity of more than a few kms -1 . 
Consequently, we decided to set the expansion velocities to zero. 
We then tried to constrain the seven remaining parameters (i, 
PA, Vrot, P, a c , fli, and EW). To reduce the number of computed 
models, we started with a qualitative estimation of the parame- 
ters from our interferometric data (especially for PA, i, a c , a/, and 
EW/) and explored the parameter space with decreasing steps to 
converge to a x 2 minimum. To check for other minima, we also 
explored the full range of possible parameters space but with 
larger steps. 

The parameter values for the best-fit model are presented in 
Table [4] The corresponding differential visibilities and phases 
are overplotted in Figs. |4]and|5] For this model, we obtained a 
reduced^ 2 of 2.0. The overall morphology and amplitude of the 
lines' profile and differential visibilities and phases are roughly 
fitted. Nevertheless, such a simple axisymmetric model cannot 
reproduce the asymmetries of the spectro-interferometric data, 
and this issue will be discussed in more detail in Section 6.3. 

6. Discussion 

6.1. The central star 

We obtain the same inclination angle i.e. i=55 + 5°, similar to the 
one determined from the continuum fit in the previous section 
and a PA. of 45 + 10°, larger than the PA obtained in the H 
and K bands (~ 20+ 5°). This is more than the value found for 
yCas (about 42°, Stee et al. 2012), and is consistent with the 
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Table 4. Parameters values for the best-fit axisymetric-kinematic 
model. 



Param. 



Value 



Remarks 



Global geometric parameters 

R* 6.5 R 

d 373 pc From von Leeuween (2007) 

i 55 ± 5° 

PA 45 ± 10° Inconsistent with Polar. 



P 



Global kinematic parameters 

450 ±50 km s" 1 See discussion 

0.5 ±0.05 Keplerian rotation 



QBrr 

EW B[y 



Bry disk geometry 

10.2+0.5 D* =1.64: 
10+0.5 A 



0.08mas 



He i disk geometry 

a H ei 7.8±1.0D* = 1.25 ±0.16 mas 

EW H d 8.+1.0A 



inference from optical Fe II emission profiles by Smith & Balona 
(2006) that the inclination should be greater than the y Cas value. 
It is also consistent with the value of 68" + 5° obtained from the 
evolutionary tracks by Ekstrom et al. (2012). 

The disk P. A. is neither perpendicular nor parallel to the po- 
larization measurement of 81° by Yudin (2001). We obtain a ro- 
tational velocity of 450 + 50 km s" 1 . Assuming that V rot , the ve- 
locity at the base of the disk is the stellar rotational velocity, we 
obtain a V sin i of 368+43 km.s" 1 . This is higher than the value 
determined by Slettebak (1982), i.e. Vsin/ 300 kms" 1 , and 
Ballereau et al. (1995), i.e. Vsin/ 318kms _1 , but compatible 
with Vsin/ = 400 + 30 kms" 1 determined by Chauville et al. 
(2001) from model fitting of the Hel 4471 and Mgll 4481 lines, 
duly corrected for the veiling effect and considering wavelength- 
dependent limb darkening inside the spectral lines. 

As already noted in Sect. 2, Zorec et al. (2005) estimate that 
the stellar mass is M/M Q = 9.6. Consequently, this gives a criti- 
cal velocity for the star of V c = 487 + 32 km s~ 1 . With this value, 
an inclination angle of 55+5°, and Chauville et al.(2005) V sini 
of 400+30kms _1 , we obtain V rot /V c =1.00+0.2. Consequently 
we can conclude that HD 110432 is rotating very close to its 
breakup velocity, like almost all interfero me trie ally observed Be 
stars. 

6.2. The disk kinematics and extension 

The disk is rotating following a Keplerian rotation, as is now 
clearly established for Be stars (Meilland et al. 2007; Meilland 
et al. 2012). Since this is the first interferometric observation 
of HD 1 10432, it is difficult to compare the geometry obtained. 
However, we can compare disk radii in the Bry with the mea- 
surements obtained in a survey of eight Be stars by Meilland et 
al. (2012). Using the same equipment and reduction algorithms 
we found that in Bry flux at least the disk size of HD 1 10432 is 
larger to those of a Col (5.5 ±0.3 R*) and a Ara (5.8 + 0.5 R*). 
Following Meilland et al. (2012), we note that a Col is a B7IV 
star with a T eff of 12963 + 203 K and a Ara a B3IV star with 
a T e ff of 18044 + 312 K. Since both are considerably cooler 
than HD 1 10432, conditions other than only the stellar effective 
temperature do influence the disk size. 

We recall that the total emitted flux in the Bry line has two 
components. One of them is due to the photospheric absorption 
of the star that underlies the circumstellar disk or envelope. The 
other component corresponds to the emission produced in the 
circumstellar environment. 



For this simple model, we can write the expression of the 
modulus of the visibility in the Bry line as 



V) 



Bry 



V+BryF *Bry + V^ nv Bry^envBry 



(1) 



where V*Bry and F *Bry represent the visibility and the flux of 
the photospheric Bry absorption respectively, while V e nvBry and 
F enV Bry are the visibility and the flux emitted by the envelope in 
the Bry line. The F tot flux is simply 



'*Bry 



+ F t 



envBry- 



(2) 



The quantity of interest in our analysis of the emitting region 
(VenvBry) can then be written by using EqQ]and Eqf2] 



Bry 



V, 



* Bi tt: 



envBry 



l - 



(3) 



Thus, changing the effective temperature will change the ratio 
^Y^- m a nontrivial way since it may increase F*Bry but also 
change the value of F enV Bry This is because modifying T e ff will 
also change the source function of the disk emitting region. The 
same argument applies to the formation of He I line and nearby 
continuum emission. 

From this paper and the cited work, the disk size does not seem 
to be correlated with the stellar spectral type. We also measured 
the disk size within the Hel emission line at 2.05 fim, which 
seems to be smaller compared to the Bry line, i.e. 7.8 vs 10.2 
D*. To our knowledge, the only other Be star with a disk mea- 
surement in the He I 2.05 fim emission line was the binary Be 
star S Sco with a disk extension of 4.5 + 0.5 D+ (Meilland et al. 
201 1), thus smaller than our measurement by a factor 1.7. In the 
Bry, Meilland et al. (2011) find a FWHM of 5.5 + ID*, again 
1.8 times smaller than the one we measured for HD1 10432. The 
difference may originate from the youth of S Scorpii disk, which 
started to form in 2000 and which was still growing in 2010 ac- 
cording to the same authors. Note that Meilland et al. (2011) 
were able to measure the disk of 6 Sco in the Bry line out to 
the inferred disk truncation radius, i.e. ~ 2.24 mas or 4.2 R* at 
periastron. But this truncation radius is impossible to determine 
for HD 1 10432 since the mass and the separation of the system 
components (assuming the star is a binary) are unknown. 
In a recent paper, Touhami et al. (201 1) have estimated the near- 
infrared continuum emission from the circumstellar gas disk of 
Be stars using a radiative transfer code for a parametrized ver- 
sion of the viscous decretion disk model. They were able to pre- 
dict the half-maximum emission radius along the major axis of 
the projected disk in the H and K bands. For HD 110432 they 
obtain a color excess E*(V* - 18/im) of 2.5 which corresponds 
to a disk density of ~ 3.1 10~ n g cm" 3 from their Fig. 1. Using 
their Fig. 2 and Table 1, it translates into a disk HWHM radius of 
~ 1.4 1-1.5 8 R* so smaller than our ~ 5.8 R*. To obtain a similar 
disk size, it seems from their Table, that the disk density is cer- 
tainly higher, i.e. ~ 1.0 10~ n g cm" 3 . On the other hand, many 
effects may explain this discrepancy: 



A larger disk inclination and a density of ~ 8 . 1 1 1 1 g cm 3 . 
The estimated radius by Touhami et al. (201 1) becomes 3.13 
R*. 

The contamination of the interferometric visibilities from the 
incoherent light of a possible companion close enough to in- 
fluence our signal at the time of the VLTI observations. In 
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that case, a correction of the visibilities from this contami- 
nation is necessary before solving for the angular diameter 
of the disk. That would considerably lower the fitting value 
of angular size. Although Mason et al. (1997) did not detect 
a companion through speckle, neither do we, and the system 
is suspected to be a high mass X-ray binary. 

• Variability: our interferometric measurements were taken at 
the beginning of 2010, whereas the AKARI IR excesses orig- 
inated before 2010. 

• The value for the density exponent adopted in the model 
(n=3) by Touhami et al. (2011) . This assumption could be 
argued to justify the higher disk size. In fact, the IR excess 
and the disk size are highly sensitive to this parameter. 

• The HWHM used in their model is the image disk radius, 
which could be slightly different from the disk radius ob- 
tained from a Gaussian elliptical disk. 

• A higher density than what the model predicts, which may 
be a clue to a new outburst episode. This could be verified 
by looking for variabilities in Ha EW, but unfortunately we 
were not able to find any Ha spectra before 2010. 

6.3. Disk inhomogeneities 

We see from Fig|4]and Fig[5]that the differential phases for these 
two lines are, first, represented well by our very simple kinemat- 
ical model, and second that they exhibit the typical "S" shape of 
a rotating disk. The amplitude of this "S" shape is ~ 10 - 15°, 
which is very similar to the differential phases obtained by 
Carciofi et al. (2009) for the Be star £ Tau, and for 6 Sco by 
Meilland et al. (2011). Our differential phases exhibit a very 
asymmetrical "S" shape with a smaller amplitude in the blue 
part of the "S" curve with respect to the central line wavelength. 
This was also the case for £ Tau. This was clear evidence for 
a one-armed spiral structure in the £ Tau disk, which may also 
be the case for HD 110432. Moreover, our visibilities across 
the Bry line are asymmetrical, with a red wing of the visibility 
systematically smaller than the blue wing. The line profiles 
are also clearly asymmetrical, with Bry and He I line profiles 
showing a V/R < 1 even if the S/N and spectral resolution for 
He I are not high enough to clearly believe this line asymmetry 
(See Fig. |6). Our simple symmetrical "toy" model is not able 
to reproduce the V/R of the observed profiles as seen from Fig. |6] 

Thus, HD 1 10432 circumstellar disk is globally rotating, with a 
larger volume and brighter region responsible for the red part 
of the line, i.e. flowing away from us, and a smaller (more com- 
pact), fainter, or more absorbed region in the blue part of the line, 
i.e. rotating in our direction. These emitting regions are also re- 
sponsible for the Bry line profiles with V/R < 1 . Moreover, since 
the differential phase is asymmetrical, it means that the photo- 
center of the emitting regions is asymmetric with respect to the 
central star (or the rotational axis). 

To account for these asymmetries, we tested two different types 
of asymmetric models. In the first one the line emission is mod- 
eled by an asymmetric gaussian disk : 

I,(r, 6) = with (r{6) = o- () (l + A cos(# - O )) (4) 

In this model the Gaussian width, given by <x, depends on the 
azimuth 6. The reference for the azimuth (0=0°) is the projected 
polar axis of the star. The parameters A and 0q set the amplitude 
and orientation of the asymmetry, respectively. 

In the second type of asymmetric model, we try to roughly 
reproduce the one-arm oscillation model presented in Okazaki 



(1997). Thus, we consider an antisymmetric spiral-like intensity 
distribution given by 

I,(r,6) = e (_ ^ ) (l + A cos(9 - 6 - nr/T)) (5) 

where A is the amplitude of the spiral perturbation of the 
Gaussian distribution, and T the spiral period. 

Examples of both intensity distributions are shown in Figure [7] 
As for the axi-symmetric model, the intensity distribution is 
then rotated on the skyplane and flattened, taking the inclination 
angle of the object into account. The use of these models 
significantly improve the fit of the data. The x 2 r obtained for the 
best-fit asymmetric Gaussian disk and the antisymmetric spiral 
are 1 .4 and 1 .6, respectively, compared to 2 for the symmetrical 
model (see section 5). For both types of models, the best-fit 
model was obtained without modifying the stellar and kinematic 
parameters from Tableland with only minor changes in terms 
of the FWHM of the intensity distribution and/or equivalent 
width of the emission line. The corresponding parameters values 
for these models, namely A, 6q, and T, are presented in Table [5] 

The two additional parameters for the asymmetric Gaussian 
disk, i.e. A and Qq are well constrained. The asymmetry is quite 
large, i.e. the largest disk extension is about 1.7 times the small- 
est one and the over intensity is orientated close to the projected 
minor axis, and mainly to the south-east. In the case of the an- 
tisymmetric spiral model, we could not obtain asymmetric visi- 
bility variations for a period of the spiral arms of less than T=20 
R*. Actually, for shorter periods, the asymmetries introduced in 
the model by the spiral averaged themselves out, considering that 
the disk FWHM in the emission lines are on the order of 4-5R*. 
Nevertheless, for longer periods, we were still able to constrain 
the two other parameters for this model. As for the asymmetric 
Gaussian disk model, we obtain an asymmetry that is roughly 
oriented in the south-east. 

Finally, the better agreement of the asymmetric Gaussian disk 
with the interferometric data is mostly because this model is able 
to produce stronger asymmetries in the visibility with almost 
symmetric S-shaped phases. Moreover, this model also repro- 
duces the visibility and phases observed along the polar orienta- 
tion better. Nevertheless, the quantity and quality of our data are 
not sufficient to fully constrain the intensity distribution in the 
emission lines. 

Regarding the possibility of a companion, we have no evidence 
of binarity from our interferometric measurements, neither from 
a modulation of the visibility modulus as a function of time or 
spatial frequency, nor from a phase modulation of the fringes. 
The question of an unseen companion orbiting HD 110432 be- 
comes an issue when studying the origin of its hard X-ray flux, 
if indeed the X-rays are formed somehow from infall in the form 
of gas stream or blobs originating in the Be star. 

7. Conclusion 

Using standard long baseline optical interferometric techniques, 
we have reported on the size and orientation of the disk of the 
first "y Cas X-ray analog" system other than the prototype y Cas 
itself. We obtained a disk FWHM of 10.1 D* in the Bry and 7.7 
D* in the He I lines. 

We have clearly detected an asymmetry in the disk in both the 
Bry and He I lines which is to our knowledge the first time that 
an asymmetry is detected on such a small spatial scale within 
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Fig. 6. Bry (Right) and He I lines from our VLTI/ AMBER observations (red line) with our symmetrical line profile fit (blue line) 
and our asymmetrical model (green line). 



Asymetric Gaussian Disk 
(Model 1) 



Antisymetric Spiral Disk 
(Model 2) 
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Table 5. Parameters values for the best-fit asymmetric kine- 
matic model. 



Parameters 


Model 1 


Model 2 




Bry disk geometry 




a Biy 


9.8+0.2D* 


10.2±0.2D* 


EW Bry 


9.5±0.5A 


10+0.5A 




He i disk geometry 






7.4+0.5 D* 


7.6+0.5 D* 


EWm 


8. ±1.0 A 


8. ±1.0 A 




Asymmetry 




A 


0.25±0.05 


0.55±0.1 




165±5° 


225 ±10° 


T 




>20 



1 



Fig. 7. Examples of pole-on (Upper figures) and projected onto 
the sky (Lower figures) intensity distributions for the two asym- 
metric classes of models tested in this paper. The central white 
circle represents the position of the central star. 



tospheres that lie between 3 R+ and 7 R* (Getman et al. 2008) 
and probably smaller than 1 R*, in order to have undetectable 
effects on our interferometric measurements. 
The disk's major-axis position angle (PA) of HD 1 10432 is 45° 
from the kinematics fit or 20° from the continuum disk emission 
fit, but in both cases is neither aligned nor perpendicular to 
the polarization measurements. Finally the disk size seems to 
be independent of the stellar parameters, as already found for 
classical Be stars by Meilland et al. (2012). 

We concluded that HD 1 10432 is a near critical rotator charac- 
terized by V rot /V c =1.00+0.2. 



the He I emission line. As a whole, the viscous disk model 
seems to be a possible scheme to match the inhomogeneities 
detected in HD 1 10432 disk since it produces non symmetrical, 
spectrally resolved visibilities and phases, as measured by 
Carciofi et al. (2009). Nevertheless, the pursuit of possible disk 
inhomogeneities requires a more dedicated set of observations, 
as done by Carciofi et al. (2009) for ( Tau. 

The disk is in Keplerian rotation and is in contact with the central 
star, as for y Cas. If HD 1 10432 has a magnetosphere, it must be 
very small, i.e. smaller than the size of classical T-Tauri magne- 
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